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^ Abstract 

r-| We consider the sensitivity of fixed-target neutrino experiments at the luminosity frontier 

^ to light stable states, such as those present in models of MeV-scale dark matter. To ensure 

Q-i the correct thermal relic abundance, such states must annihilate via light mediators, which in 

(— I turn provide an access portal for direct production in colliders or fixed targets. Indeed, this 

'~~' framework endows the neutrino beams produced at fixed-target facilities with a companion 

^ 'dark matter beam', which may be detected via an excess of elastic scattering events off 

^ electrons or nuclei in the (near-) detector. We study the high luminosity proton fixed-target 

QQ experiments at LSND and MiniBooNE, and determine that the ensuing sensitivity to light 

^O dark matter generally surpasses that of other direct probes. For scenarios with a kinetically- 

mixed U(l)' vector mediator of mass my^ we find that a large volume of parameter space is 

Q excluded for mDM ~ 1 — 5 MeV, covering vector masses 2mDM ^ f^v ^ f^-q and a range of 

^-H kinetic mixing parameters reaching as low as /t ~ 10~^. The corresponding MeV-scale dark 

matter scenarios motivated by an explanation of the galactic 511 keV line are thus strongly 

. ^ constrained. 
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1 Introduction 

While the empirical evidence for dark matter (DM), through its gravitational effects in 
astrophysics and cosmology, derives from many sources and ranges over many distance scales, 
the search for any signature of its non-gravitational interactions remains one of the focal 
points of research in particle physics. Thermal relic weakly-interacting massive particles 
(WIMPs), predicted or otherwise introduced into many extensions of the Standard Model 
(SM), represent an appealing dark matter candidate. In particular, though the WIMP mass 
scale and couplings are only weakly constrained by the requirement of the correct relic 
abundance, the characteristic weak-scale parameters of the paradigmatic WIMP fall into a 
range that offers hope for the direct discovery of non-gravitational DM interactions in the 
laboratory. However, in recent years, considerable attention has been paid to particle physics 
models that deviate from the minimal idea of a single WIMP with weak-scale interactions 
with the SM. Possibilities of both light dark matter candidates and/or light mediator particles 
have been explored fT-17 with the motivation of tying various anomalous experimental 



signatures to the annihilation, scattering or decay of dark matter. While most anomalies 
will likely find other explanations, the expanded mass range for WIMP candidates and 
mediators opens a number of new experimental avenues, which go beyond the characteristic 
direct detection strategies for a minimal (weak-scale) WIMP. 

In this paper, we revisit a class of MeV-scale dark matter models, originally designed 
to explain the unusual strength and morphology of the 511 keV emission observed from the 
galactic center with annihilating dark matter [ll. MeV-scale models of thermal relic dark 
matter require the existence of a light mediator [2]-[9] , so that cosmological freeze-out occurs 
as a rescaled version of conventional WIMP freeze-out with reduced mass and temperature 
scales. The existence of light mediators, and thus a more complex light hidden (or dark) 
sector, naturally stimulates interest in the low and intermediate energy particle physics man- 
ifestations (2, 10 14 . Certain classes of flavor-conserving light states with lifetimes below 1 



second are often immune to a variety of astrophysical, cosmological, and collider tests, even if 
their interactions are larger than the characteristic weak rate. Specifically, in scenarios where 
states in the hidden sector are in the hadronic mass range, and have a lifetime longer than 
other hadronic states which undergo weak decays, fixed target experiments with detectors 
10-lOOOm from the target, as in modern long-baseline neutrino experiments, can provide 



complementary sensitivity to colliders 15 17 . A rather striking consequence of models with 



light (sub-GeV) dark matter is the production of a high intensity 'dark matter beam', gener- 
ated as dark matter particles are pair-produced as a result of the proton-target interactions 
and boosted along the proton beam direction [16j. The scattering of light dark matter in the 
(near-)detector would then generate an additional source of neutral-current-type scattering 
events (see, e.g. |!l8|). This prediction implies that a direct search for MeV-scale stable dark 
matter is possible at experiments at the luminosity frontier, which is the focus of the present 
paper. 

To motivate the importance of fixed target facilities in this low mass regime, we recall p6] 
that given such a hidden sector, assumed neutral with respect to the SM gauge group, we 



can parametrize the interactions as follows, 

(nidi) (nidi) 
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where O denotes SM and new physics (NP) operators of canonical dimensions di and d2, 
and A is a cutoff scale presumably at a TeV or above. Light, long-lived, hidden sector 
states can be studied at high- luminosity fixed target experiments, where the production 
cross-section mediated by an interaction ([I]) of dimension di + d2 = 4: + n (with n > 0) 
typically scales as a ~ £'^"~^/A^"'. Inserting typical numbers for the attainable luminosities 
and typical energies at high-energy colliders and proton fixed-target machines respectively, 
leads to an interesting comparison in the total production count (denoted A^) of neutral 



GeV-scale states 16 
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It is apparent that for a marginal interaction, n = 0, the production rates at fixed targets may 
be sufficiently advantageous to easily counteract the low geometric acceptance of a detector 
placed some distance from the target, and such facilities can provide the dominant level of 
experimental sensitivity. The set of relevant or marginal interactions forms a small, but 
generically the most important, subset of interactions in ([I]) known as SM portals 19 22]: 



0-^' = F^^,H^H,LH, (3) 

where F^, H and L are the hypercharge field strength, and the Higgs and lepton doublets. 
The operators (IS]), denoting respectively the vector, Higgs, and neutrino portals, allow a 
coupling of the SM to (SM neutral) new physics at the renormalizable level. 

This paper aims to explore the sensitivity of existing (and future) experimental infras- 
tructure for long-baseline neutrino experiments to light dark matter which forms part of 
a hidden sector interacting with the SM through the portals ^. At fixed targets, these 
interactions entail the production (along with the neutrinos) of a boosted dark matter beam 
through the generation and subsequent decay of GeV-scale mediators. As we will show below, 
existing data from high-luminosity experiments such as LSND and MiniBooNE already im- 
poses stringent constraints on viable scenarios of MeV-scale dark matter, due to the limits on 
neutral-current-like scattering events off electrons and nuclei in the detector. While the idea 



of searching for exotics using fixed target facilities is certainly not new (see e.g. 18 23 26]), 
long-baseline neutrino facilities introduce a particular advantage for probing stable states in 
that the large mass of the (near-) detector can be utilized to observe scattering rather than 
just the results of a decay. We will argue that these facilities already provide the dominant 
constraints on many models of this type, and specifically those utilizing the vector portal 
C^SM = F^i, which are generally less constrained in other ways. In particular, for scenarios 
of MeV-scale dark matter which aim to explain the galactic 511 keV line, these constraints 
are generically more stringent than those derived from rare meson decays [2]. 

The rest of this paper is organized as follows. In section [2| we describe and motivate one 
of the most viable classes of MeV-scale hidden-sector dark matter, which interacts with the 
SM via the vector portal, and the parameter constraints from astrophysical and cosmological 



data. In section |3} we explain how the parameter space of this model may be probed at fixed 
target neutrino oscillation experiments and follow this with an analysis of the sensitivity of 
the LSND and MiniBooNE experiments. We conclude by contrasting this sensitivity with 
other limits that can be placed on MeV scale states, and explore possibilities for future 
progress, in section |4j 

2 Light MeV-scale thermal reUc dark matter 

The viability of thermal relic dark matter with a mass in the MeV-GeV range, well below 
the Lee- Weinberg bound, seemingly rests on the presence of a light hidden sector with states 
which can mediate annihilation (2l[3l|6]. Moreover, various phenomenological constraints |6 
suggest that the most viable scenarios are those in which the hidden sector is uncharged 
under Standard Model symmetries. This naturally leads us to the portal interactions ([s]), as 
the primary means of probing these sectors at low energies. 

For a thermal relic dark matter (TRDM) candidate in the MeV mass range, the dominant 
decay channels will lead to e^e~, with direct annihilation to photons and neutrinos often 
suppressed (a reduced coupling to neutrinos being a phenomenological constraint to eliminate 
drastic softening of the supernova neutrino spectra). The fact that annihilation of light 
TRDM generically produces positrons naturally implies that galactic observations are a 
significant source of constraints. Indeed, MeV-scale models were initially motivated by the 
511 keV line observed from the galactic centre [I], and recently mapped out in considerable 



detail by INTEGRAL/SPI [27,28 . However, it is important to emphasize that, independent 
of any attempt to explain its source, the magnitude of this fiux provides quite a significant 
constraint on light dark matter models in this class. To see this, we parametrize the fiux $ 
observed by INTEGRAL/SPI in the form. 
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where A^e+ is the number of positrons per decay and ^DM/^m ~ (lpbn)/(crf)fo. The result 
depends crucially on the annihilation rates at freeze-out (erf )& and in the galactic centre 
{av)gai- We see that if annihilation is dominantly to e^e~, any MeV-scale dark matter 
candidate with relic abundance close to Qm should have an annihilation rate in the galaxy 
suppressed by several orders of magnitude relative to the rate at freeze-out. This is quite a 
strong constraint, and as emphasized in [3|[6] tends to single out one class of models as the 
most viable. We summarize the issues below for the case of interactions mediated via the 
vector and Higgs portals, with a light hidden sector containing a dark matter state and a 
U(l) mediator V or singlet scalar mediator S respectively (gI: 



(/) Vector portal, /tidm < "mv'- A scalar DM candidate has p-wave annihilation which 
satisfies Q since v ~ 10~^ in the galaxy, and thus is viable for sub-percent mixing via 
the portal coupling. 

(X) Vector portal, ttidm > "^y: This implies s-wave annihilation, and thus the 511 keV fiux 
limit Q can only be satisfied with a highly subdominant component of dark matter. 



(X) Higgs portal, tudm < f^s'- Annihilation is suppressed, and would require 0(1) mixing 
via the Higgs portal which is ruled out for example by K and B decays. 

(X?) Higgs portal, r/iDM > "^5: A fermionic DM candidate has p-wave annihilation which 
can satisfy (ll]), but needs a high degree of tuning to avoid limits on K decays with 
missing energy. 



This phenomenological analysis, described in more detail in l6|, motivates models inter- 
acting via the vector portal as the most natural setting for light MeV-scale TRDM, with the 
dark matter state being the lightest in the hidden sector. In the next subsection, we will 
outline this model in more detail and then move on to explore how it may be probed in fixed 
target experiments. 

2.1 Light dark matter with a U(l) mediator 

We consider a hidden sector, charged under a U(l)' gauge group, with a vector portal coupling 
to the SM via kinetic mixing (see e.g. |6]). We also assume that the U(l)' is spontaneously 
broken at a low scale by a Higgs' sector, leading to a mass for the vector mediator V^. The 
relevant low energy Lagrangian takes the form, 

^v,x = -\y'u + \ri\yl + >^y.Q^F^'' + i^^xi' - ^\\xf + A', (s) 

where x is the complex scalar dark matter candidate, taken to be stable due to a suitable 
Z2 parity, and the U(l)' covariant derivative is D^ = 9^ + ze'V"^. All of the kinetic terms and 
interactions involving the Higgs' are included in the £/j/ term and will not play a role here, 
although we require that the full scalar potential is such that the physical Higgs' is more 
massive than x- In general, V mixes kinetically with the hypercharge gauge boson, but at 
low energies we can ignore the induced coupling to the Z, and in ([s]) we have rescaled the 
kinetic mixing parameter k so that it just reflects mixing with the photon. 

The model contains four parameters; the masses m^ and my of the dark matter candidate 
and the vector mediator, the U(l)' gauge coupling e', and the kinetic mixing coefficient k. On 
requiring that x comprises the majority of dark matter, the constraint on its relic abundance 
allows us to fix one relation between these four parameters. The primary quantity here is the 
annihilation rate, which is given in general by the diagram on the left of Fig. [1} In practice, 
we will be in a regime here where the branching is predominantly to an e^e" final state. In 
the limit of small mixing, and dropping a small correction proportional to rrig, the rate for 
annihilation is given by p3]. 




(a.)_ :. 3 X 10-- cm^ x — (^^) x-— xJl--f -—^-^ . (6) 



Noting the p-wave suppression, with v ~ 0.3 at freeze-out, we observe that the WMAP con- 
straint on the relic density f^DM^^ ~ 0.1 ~ (0.1pbn)/(o"f)fo imposes the following restriction 
on the model parameters, 
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Figure 1: Tree-level annihilation (left) and scattering (right) of scalar dark matter in the U(l)' hidden 
sector. 

where to simplify the presentation we have taken ml <^ rn^ <^ m^, relations which are satis- 
fied up to 0(25%) in the parameter regimes studied here. However, we use the more precise 
constraints from ^ in the subsequent numerics, reducing the number of free parameters 
to three, which we will take to be {m-^, my, k}. Note also that the p-wave suppression of 
annihilation for low velocities allows this process to satisfy the 511 keV fiux constraint (|4]), 
as alluded to above, as well as the CMB constraints on dark matter annihilation 1291. 



A rotation of the annihilation diagram describes scattering off electrons and nucleons as 
shown on the right of Fig. [T| and provides a means for detecting the presence of light DM 
in the galactic halo, at least in principle. However, in practice MeV-scale dark matter only 
has a characteristic kinetic energy of (9(eV) in the Earth's rest frame, leading to a recoil 
in nuclear scattering which is well below the detection threshold for the current generation 
of underground direct detection experiments. However, this problem could be circumvented 
if dark matter were first boosted to w ~ 1 before it undergoes scattering. Indeed, for 
sufficiently low masses, it is feasible to produce a dark matter beam in collider or fixed target 
experiments that may see MeV dark matter in the ultrarelativstic regime, E^ 3> m^. In 
particular, modern long-baseline neutrino facilities are ideal for this purpose, having a high 
luminosity, and also large (near-)detectors, which can be used to search for dark matter 
elastic scattering. We will turn to this possibility, and the sensitivity that can be attained, 
in the next section. 



3 Production and scattering of the dark matter beam 

Within the hidden sector scenario outlined in the previous section, and for sufficiently small 
mvi the following chain of processes can produce a dark matter beam at a fixed target 
experiment: 

1. p -|-p — > X + 7r°,r7 

2. 7r°,?7 ^7 + 1^ 

3. V -> 2x 

Depending on the beam energy and form of the target, the relevant decay lengths ensure 
that this entire sequence of events will occur either inside the target itself or in the subsequent 
decay volume. Depending on my, the dominant production mode will be 7r° or r} decays. 




and we have focused on this subset of hadronic states due to their large branching fraction to 
photons. While vr'^'s dominate production for mv < "^vr, the addition of the r] mode allows 
access to a larger range in mv, and consequently m^. In both cases, the branching ratio 
to V is proportional to that of the radiative decays of the mesons to two photons, though 
suppressed by k^ and phase space factors related to the ratio of mv to m^ where = tt, r/, 

Br^^^y ~ 2fi;^ 1 ^ i dl^^^^. 

For the case of vr^ decays Br^o_j.^^ ~ 1, while for r] decays Br^_j.^^ ~ 0.39. 

Given that we require k <^ 1, it follows that V generically decays within the hidden sector, 
Btv^2x — I5 ^^d the ensuing dark matter beam then propagates along with the neutrino 
beam. For the range of k values considered here, it has a weak-scale scattering cross-section 
with normal matter and may be detected through neutral current-like processes, either with 
electrons e -|- x ~^ e + x, or nucleons, N + x ^^ ^ ~^ X- ^^ order to probe this scenario, we 
will utilize the results of LSND and MiniBooNE, which have two of the largest datasets and 
importantly have published analyses on neutrino elastic scattering, which DM scattering will 
closely mimic. Note that due to their respective beam energies, both LSND and MiniBooNE 
are sensitive to electron scattering, while only MiniBooNE is sensitive to elastic scattering 
off nucleons. 

3.1 Dark matter beams at LSND 

We now probe the parameter space of the model by calculating the number of dark matter 
neutral current-like elastic scattering events that would be expected at the LSND experiment, 
-^events, and comparc it to the total number of elastic (neutral and charged current) scattering 



events off electrons actually observed 30 . At LSND, pions were produced by impacting an 
800 MeV proton beam onto either a water or high-Z metal target I31j. The LSND experiment 
provides the largest fixed-target sample of pions currently available, and has the potential 
to provide the most stringent limits on the model parameter space for the range of my's to 
which it is sensitive. 

The overall normalization of the event rate at LSND is dictated by N^o, the total number 
of neutral pions produced over the lifetime of the experiment. In practice, we can approx- 
imate iV^o by equating the 7r° production rate with that of tt"*", on the grounds that the 
measured 7r° and 7r+ production rates in proton-nucleon collisions differ by 0{1) factors (see 
eg. [32]). We estimate Nt^+ by working backwards from the neutrino flux reported by the 
collaboration, and as the majority of the neutrinos were products of 7r+ decays at rest. 

Here $;y = 1.3 x 10^^ u cm~^ is the neutrino flux over the lifetime of the experiment, 
y4det — 2.5 X 10^ cm^ is the area of the detector facing the target, and {dQiah/4:7!')u ~ 3 x 10~^ 
is the fraction of the solid angle subtended by the detector relative to the target. 



We employed a Monte Carlo simulation to determine the dark matter flux incident on 
the LSND detector. Pions were generated in the momentum ranges expected by LSND 
over an array of possible angles. According to the appropriate branching fractions, the 
subsequent decays to 7r° — )■ V'j and V — )■ xx^'s, were simulated and the trajectories of the 
x's were then checked to determine if they intersected with the detector. A re-weighting 
technique was then used to weight each trajectory according to the momentum and angular 
distribution of the initial 7r°. We assumed that this distribution was similar to the production 
distribution of vr+'s and used the parameterization of the production cross-section by Burman 
and Smith |33J. It was also necessary to account for the fact that the pion production 
distribution was not constant throughout the lifetime of the experiment, as LSND made 
use of two different targets (3l], and the resulting normalized (and azimuthally symmetric) 
distribution, which we denote f^^{0,p), is a weighted average. 

With the pion distribution in hand, the simulation determined the number of x's which 
reach the detector, along with their energies and the distance travelled through the detector. 
To determine the expected number of elastic scattering events, we modelled the detector as 



a cylindrical tank filled with mineral oil CH2 (see 34 for further details). The scattering 



cross section for ex -^ ex, assuming E ^ 771^, takes the form 

dEf a ^ E^{ml + 2m^EfY ' ^ ' 

where E is the energy of the incoming dark matter particle and Ej is the energy of the 

I I / 2 

scattered electron. We use (7) to replace the ratio ^-^ with a function of my and m^. The 
number of elastic scattering events of dark matter off electrons can now be schematically 
represented as follows, 

A^e™X = ^^ ^ ^-" ^ Br,0^7V X eeff X J] Liae^^ex(^i)/r(^*,Pi)A», (H) 

i 

where Ue — 5.1 x 10^^ electrons/cm" is the number density of electrons in mineral oil, 
Br^o^^y is calculated in (jsl), and e^s — 0.19 is the electron detection efficiency within LSND's 



elastic scattering analysis 30 , that we will use for comparison. In order to map out the 
full production distribution, the remaining sum is over all simulated trajectories, where L 
is the distance travelled through the detector and aex-^ex is the integrated cross section 



from (10) for outgoing electrons with recoil energies between 18 MeV and 50 MeV. These 



cuts are again chosen to match those of LSND's ve elastic scattering analysis 30 . Finally 
A = 6pT,g696(j)/{27i) reflects the step size in solid angle and pion momentum, corresponding 
to each simulated trajectory. 

The results of the simulation, for a range of values of my up to the pion threshold are 
shown in Fig. [21 where we have plotted the expected number of events (from ( 11 )) that would 



be detected by LSND in the K-my parameter space for dark matter with m^ = 1 MeV. The 
sensitivity is consistent with the earlier analysis in [16|. In addition, the solid black line 
delineates the strong coupling boundary, below which a' > 47r. While this marks the regime 
where our perturbative calculations cease to be reliable, constraints on self-interaction can 
be somewhat stronger than this. LSND only observed (9(300) beam-on events passing their 
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Figure 2: Expected number of elastic scattering events of dark matter off electrons at the LSND detector 
for m^ = 1 MeV. The regions show greater than 10 (light), 1000 (medimn) and 10^ (dark) expected events. 
The area below the black line corresponds to a' > An. 



cuts, of which around 200 were expected to be elastic scattering events due to neutrinos 30 . 



Thus, in hght of the fact that DM scattering events should pass the same cuts, we can easily 
exclude the medium and dark regions of Fig. [2| for which LSND should have observed far in 
excess of a few hundred events. These limits could be further improved by a spectral analysis 
of the recoiling electrons, as the energy of the dark matter beam is considerably higher than 
that expected for the neutrinos from pion and muon decays at rest, and correspondingly 
recoil electrons from dark matter scattering would be more energetic than those from the 
neutrino beam. With such an analysis, one could plausibly compare the expected number of 
events with those due to the far rarer neutrinos produced by pion decays in flight, of which 
LSND observed (9(10). However, even without the extra sensitivity that spectral information 
would provide, we can exclude the majority of the parameter space for 10 MeV < mv < m.^ 
(for m^ ~ 1 MeV) through a combination of the strong coupling condition and the expected 
number of NC-like scattering events. We will next look at the sensitivity of the MiniBooNE 
experiment, which can probe parts of the model parameter space with larger m^ and my. 



3.2 Dark matter beams at MiniBooNE 

The MiniBooNE experiment made use of an 8.9 GeV proton beam impacting a Be target [351, 
and thus can produce a more energetic dark matter beam. The process for calculating the 
number of neutral current-like elastic scattering events is quite similar to that described for 



LSND in section pA] though with some differences as outlined below. 

As for LSND, the number of vr^'s produced, A^^o, was approximated by the number of 



charged pions produced over the hfetime of the experiment, A^^. As most of these pions 
do not decay at rest, we account for the forward boost by replacing (ifiiab in pi) with the 
fractional solid angle in the pion centre of mass frame, dQcm — I'^d^iah for small angles. The 
average vr"*" energy was 1.12 GeV [35], which corresponds to 7 ^ 8. Using A^et — 1-2 x 10^ 
cm^, $^ = 3.35 X lO^^i^cm-^ and (dni^^/ATr)^ = 3.2 x 10-^ leads to N^+ ~ 1.6 x lO^". 
However, this is a significant over-estimate, due to the influence of a magnetic focusing 
horn, used to allow the experiment to run in either neutrino or anti-neutrino mode. The 
influence of the horn is easy to isolate, though, as the beam was run for a short period with 
the horn turned off. The neutrino flux was observed to drop by a factor of six during this 



period 36 , and we decrease the value of A^^+ accordingly, arriving at the following estimate 
for 7r° production at MiniBooNE: A^^o ^ 2.6 x 10^^. 

It is apparent that the number of neutral mesons produced at MiniBooNE is nearly three 
orders of magnitude lower than the number at LSND, but this is compensated to a significant 
extent by the large forward boost that the beam acquires, which tends to enhance the number 
of x's whose trajectories intersect the detector. Moreover, the MiniBooNE proton beam is 
of considerably higher energy than that of LSND, and as such is capable of producing ?7's 
in significant quantities. With a mass of 547.8 MeV, including rj production allows us to 
greatly extend the range of V masses that can be probed using these fixed target neutrino 
experiments. In addition, larger values of m^ are also accessible, and we have calculated 
MiniBooNE's sensitivity to m^ = 50 MeV in addition to m^ = 1 MeV. To estimate the r] 



production rate, we make use of some early experimental data 32 , which indicates that in 
the appropriate energy range. 

We use this ratio to normalize the number of r^'s produced over the lifetime of the experiment 
to A^ttO, and arrive at an estimate of A'^^ ~ 9 x 10^'''. In order to determine the sensitivity of 
MiniBooNE to the model, we will combine the results for dark matter from both 77 and 7r° 
decays. 

The Monte Carlo simulation follows similar lines to that described for LSND. The nor- 
malized distributions for 7r° and rj production in this case were approximated by averaging 



the Sanford and Wang fits for vr"*" and vr^ production used by MiniBooNE 35 , which we 
denote f^o^{0,p) ~ f^^{0,p)- The number of expected electron scattering events then takes 
the schematic form, 

N^l,^ = n,x J2 iV^Br^^^v'X^L,a,^^ex(^«)/r(^-P*)^- (13) 



where we have used the same notation as (11), and for MiniBooNE n^ = 5.1 x 10^^ 
electrons/cm^. A^^o and A^^ are given above and provide the overall normalization, while 



the electron scattering cross section is given in (11). Ultimately, we find that the ability to 
probe the model through neutral current-like elastic scattering events with electrons is still 
somewhat weaker than LSND for mv below the pion threshold, but crucially it can extend 
the sensitivity range for mv up to ?7i^. We should note that there is currently no published 
experimental analysis for elastic scattering with electrons at MiniBooNE, so while we will 
determine the potential sensitivity of the experiment, this will be overly optimistic as there 
are no cuts imposed on the electron recoil and we have ignored other efficiency factors. 
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Figure 3: Expected number of neutral current-like dark matter electron scattering events at the MiniBooNE 
detector for ruy, — 1 MeV. The regions show greater than 10 (light), 1000 (medium) and 10^ (dark) expected 
events. The plot on the left shows dark matter resulting from 7r° decays, while the plot on the right combines 
dark matter from both tt" and 77 decays. The area below the black line corresponds to a' > An. 

In order to provide a direct comparison between the sensitivity of MiniBooNE and LSND, 
we first present the estimated number of neutral current-hke elastic dark matter electron 
scattering events for m^ = 1 MeV dark matter produced via pion decays in the left panel 
of Fig. |3} The shape of this plot is very similar to that found in Fig. [2| but MiniBooNE's 
sensitivity is down by an order of magnitude, even in the absence of cuts. However, in the 
right panel we present a similar plot incorporating the contribution from r] decays. The bump 
in A^events at low my represents the pion contribution, and while it drops by about an order 
of magnitude for higher masses, the plot clearly illustrates the utility of the MiniBooNE 
dataset in providing sensitivity all the way up to my ~ 0.5 GeV. 

While there is currently no published analysis of electron elastic scattering, MiniBooNE 
has recently published a full analysis of NCE scattering off nucleons 



36 



which may lead to 

a substantial boost in sensitivity relative to electron scattering for certain mass regimes. In 
order to utilize this channel, we require the differential cross-section for neutral current-like 
elastic scattering between dark matter and nucleons (see Fig. [I]). This process is somewhat 
similar to the (vector part of) Z-mediated neutrino-nucleon elastic scattering (see e.g. [37]), 
and we obtain the following differential cross section. 



da 



XN^XN 



a'K' 4vra2 [Fl^{Q^)A{E, E^) - \Fl^{Q')B{E, E, 
X 



a 



{m\ + 2mN{E - E^y {E^ 



(14) 



mt 



where E and E^ are the energies of the incident and outgoing dark matter particles, re- 
spectively and Q^ = 2ini^{E — E^) is the momentum transfer. The functions A and B are 
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defined as: 



A{E,E^) 
B{E,E^) 



2mNEE^- m\{E - E^), 

{E^ - E) [{E^ + Ef + 2mN{E^ - E) - Ami] 



(15) 
(16) 



The cross section holds for both neutrons and protons so long as the appropriate nuclear 
form factors are used. We make use of the following monopole and dipole form factors 



'l,Ar 



qn 



(1 



■ QV^l) 



Fin — 



Kn 



[I + Qym 



'n) 



(17) 



1.9. Equation (14) only describes the 



where qp = 1, g„ = 0, while Kp = 1.79 and k^ 
scattering cross section for free nucleons, but in reality, dark matter will be scattering off 
nucleons bound in a carbon nucleus or one of the two hydrogen nuclei of a CH2 molecule. 
Following the MiniBooNE analysis 



36 , we write the effective differential cross section as 



follows, 



dE^ 



7 



CpfiQ 



+ ^Cp,m 



da 



XP^XP 



dE^ 



I r< /r)2\"''^xn^xn 



'li 



where the C's describe MiniBooNE's relative efficiencies for detecting scattering off one of 
the protons or neutrons bound in a carbon molecule or one of the protons making up the 
hydrogen nuclei. The efficiencies are dependent on the momentum transfer of the scattering, 
but are quite close to unity for Q^ G [0.4, 1] GeV^ 36 . We have not made a distinction 



between bound and free proton scattering cross sections in ( 18 ) or the numerical analysis 



The final expression for the expected number of NCE-like nucleon dark matter scattering 



events at MiniBooNE is very similar to (13). Approximating the chemical composition of 



the mineral oil used in the detector as pure CH2, we obtain 



yyMB 



■Nx 



14ncH2 X Ceff X ^ iV^Br^^^^y ^ L 



off 



iC^Nx^Nx 



{E,)f'^^{e,,p,)A, 



(19) 



(f)=TT^\ri 



where ncH2 = 6.4 x 10^^ molecules/cm^ is the number density of CH2 in the mineral oil used 
at MiniBooNE, while ees — 0.59 is the detection efficiency 38 for events within the specific 



36 



fiducial volume and momentum transfer cuts imposed in the MiniBooNE NCE analysis 
that we adopt here to allow for a direct comparison. The momentum transfer cut of 0.1 
- 1.6 GeV^ determines the range over which (18), which weights over proton and neutron 



scattering, is integrated to produce the effective cross-section cr^^^jvx- Note that the lower 
cut at 0.1 GeV^ means that there is no sensitivity to coherent nuclear elastic scattering, and 



our nucleon-level treatment should be reliable. The remaining factors in (19) are defined as 



in (13). 



The results for DM nucleon scattering at MiniBooNE are shown in the left panel of Fig. 
|4j In contrast to the situation for NCE neutrino scattering, the nucleon scattering sensitivity 
here is slightly weaker than for electron scattering at low V masses, though it does improve 
as we increase the V mass. This apparent suppression of the nucleon scattering cross-section 
can be understood as follows. The usual NCE-type enhancement naturally emerges here in 
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Figure 4: Expected number of neutral current-like dark matter nucleon scattering events at the MiniBooNE 
detector. The regions show greater than 10 (light), 1000 (medium) and 10^ (dark) expected events. The plot 
on the left is for m^ — 1 MeV, while the plot on the right is for rriy^ — 50 MeV. The area below the black line 
corresponds to a' > 47r, while the dashed curve indicates the total number of (background) neutrino events 
observed. 

the rriy ^ Q^ limit, where the relative nucleon vs electron scattering cross-sections scale 
as (XNj^^N^/aex^ex ~ rnN/fne ~ 10^. In contrast, the characteristic regime covered here 
is Q^ > my and this large enhancement for nucleon scattering is lost. Nonetheless, we 
see that as my increases the sensitivity increases, and actually peaks around mv ~ 300 
MeV for events due to rj decays. In practice, fully exploiting this strong underlying sensi- 
tivity at MiniBooNE would require a careful analysis of the spectral information, in order 
to isolate DM scattering events within the total elastic scattering dataset which contains 
approximately 9.5 x 10"^ events (of which roughly 65% are expected to be actual NCE scat- 
terings) [36] . While all these events are necessarily similar, and we rely on this to conclude 
that DM scattering would pass the signal cuts, the characteristic DM beam energy should 
differ from the characteristic energy of the neutrinos, so separating the signal from the neu- 
trino 'background' seems feasible at some level of precision. However, even without a more 
sophisticated analysis of this type, the fact that the strong coupling condition a' < An be- 
comes increasingly restrictive at higher V masses, allows us to rule out almost the entire 
parameter space of 1 MeV dark matter for 2m-^ < my < m^^^, particularly when these results 
are combined with the exclusion region from LSND. 

Independent of any motivation for such models as explanations of the galactic 511 keV 
line, it is interesting to explore how far the sensitivity reach extends in DM mass. It is 
in this spirit that we show the result of repeating the previous analysis for dark matter 
with rrix = 50 MeV in the right-hand plot of Fig. Ill We find that there is still significant 
sensitivity of MiniBooNE to this dark matter mass range, while the bounds imposed by the 
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strong coupling condition are somewhat weaker. Thus a significant portion of the parameter 
space is still allowed by experimental data, under the conservative assumption that the 
(9(10 — 1000) events for the parameter region above the strong coupling boundary could be 
hiding in the large NCE dataset. While not shown here, this calculation was repeated for 
rriy^ = 100 MeV leading to similar, though slightly weaker, bounds. 

4 Discussion and Conclusions 

There has been a recent resurgence of interest in MeV-to-GeV scale phenomenology, in 
which considerable attention has been devoted to the search for light mediator particles. 
If such mediators are the lightest members of an extended hidden (or dark) sector, it is 
then guaranteed that if produced they will at some point decay back to the SM. An explicit 
example is a new light U(l)' vector coupled to the SM via the vector portal. A number of 
search strategies are based on the production of such a vector with its subsequent decay to 
electron-positron pairs (see e.g. fl5]). However, these search strategies are severely limited 
if, instead, the mediator is not the lightest hidden sector state and decays predominantly 
to some other light states. In this case, the scattering of those light states in a detector 
spatially separated from the production point represents perhaps the most efficient search 
strategy. Moreover, owing to the potentially large production rate, and the existence of 
large volume detectors, proton fixed-target facilities focusing on neutrino physics appear to 
be an ideal means for exploring these scenarios. In this paper, we have given a detailed 
example of how such a search could be undertaken, with MeV-scale DM models as our 
primary focus. In turn, we have found that the most natural/economical model realization 
of TRDM in the MeV range - with vector portal mediation - is under severe pressure from 
the experimental constraints imposed by the absence of any non-standard elastic scattering 
signal at MiniBooNE and LSND. In particular, in the limited mass range m^ ~ 1 — 3 MeV 
that is required by other 7-ray limits [39] for any putative explanation of the galactic 511 keV 
line, these constraints rule out this candidate for most values of the vector mass below the 
rj threshold. 

In the remainder of this section, we will comment briefiy on other sources of constraints 
on these models, which are generally less sensitive than the fixed target probes discussed 
here, and also mention some future directions. 

• Additional particle physics constraints: To illustrate the advantage of constraints derived 
from the DM beam analysis, we can compare them with those that follow from meson decays 
and precision QED measurements. The anomalous magnetic moment of the muon is a very 
sensitive probe of sub-GeV physics (see e.g. (2|[II]). In the vector portal scenario studied 
here, additional positive contributions to g — 2 originate from a muon-mediator loop, and the 
current sensitivity corresponds to mixing angles in the range k ~ 10~^ — 10~^ [111. From the 
results of this paper (Figures 2 and 4), one can see that if the vector mediator is unstable 
with respect to decays to lighter MeV-scale states, the constraints derived from LSND and 
MiniBooNE are far more stringent, and in particular do not allow us to ascribe the current 
discrepancy between the measurement and the SM calculations of {g — 2)^ to a light vector 
mediator with k ~ O(10~^). 
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The vector isosinglet nature of the minimal vector portal provides a natural defense 
against excessive flavor-changing effects in the meson sector. Nonetheless, the two-body 
decay K^ -^ n^V followed by the invisible decay of the mediator will give an additional 
contribution to the tiny, but measured, branching fraction K~^ — )■ vr"*"^. Using the results 



of 11 , we estimate that the maximal contribution of the vector mediator to the decay rate 
cannot exceed Byk^ttV^tt^ ~ 3 x 10^^° x (k/IO^^)^. Given that this decay is observed 
with about this branching ratio, we conclude that the sensitivity to k, does not exceed 10"^, 
which is again inferior to the constraints derived in this paper. Finally, direct experimental 
constraints on the missing energy events e~^e~ — )■ jV — )■ 7^ at B-factories were discussed 
in 10 . The sensitivity of this process to n currently does not exceed 10~^, as such a search 



is very difficult to implement. 

• Astrophysical and cosmological constraints: We should also address the question of as- 
trophysical constraints on light mediators and MeV-scale DM. As we have discussed, the 
vector portal model considered here naturally produces a galactic 511 keV flux below or 
approaching the observations by INTEGRAL/SPI, but there are other possible sources of 
astrophysical or cosmological constraints. For example, in the early universe the lower end 
of the mass range for MeV DM implies that it will freeze-out quite close to the BBN epoch. 
The annihilation products can then have an effect on the light element abundances. This 



question was analyzed in 40 , with the result that this effect is rather small for rriy > 1 



MeV. Another possible source of constraints is the impact of DM annihilation on the CMB 



in the late universe 29 . Such constraints can be quite stringent for MeV-scale DM that has 
an s-wave annihilation rate, as the annihilation products have energies in an ideal range to 
heat the IGM. Nonetheless, in the present case these constraints disappear as the p-wave 
annihilation is highly suppressed in this epoch where DM is very cold. 

A potentially promising direct source of (semi-)relativistic MeV-scale WIMPs are certain 
extreme astrophysical environments. Indeed, the maximum attainable energies in the core 
of supernova explosions are in excess of 10-20 MeV. However, given that the interaction 
rate of MeV dark matter with the electrons is appreciably higher than the weak rate, as is 
the case for the parameter range studied here, any DM states produced will be thermalized 
and trapped in the core, and thus will not lead to additional constraints via new energy 
loss mechanisms. Given its coupling to the SM via the vector and/or Higgs portals, it is 
also guaranteed that the neutrino spectrum is not degraded and remains consistent with the 
energy range suggested by the detection of the SN1987a neutrino signal. To escape the core, 
DM would have to diffuse the same way neutrinos do. However, while the neutrino energy 
spectra can be characterized as almost thermal with a 'temperature' in the range of 5 to 
10 MeV, the corresponding MeV DM 'temperature' will have to be smaller due to its larger 
interaction rate with electrons and positrons. It is then guaranteed that the SN freeze-out 
'temperature' for x will be T^ < me- A population of energetic light DM particles created 
by past SN explosions is in principle detectable via their interaction with e.g. electrons and 
residual ionization. However, given that the diffuse SN neutrino background still remains 
undetected, discovering MeV scale particles this way represents a serious technical challenge. 

• Future progress: One interesting avenue for extending the sensitivity reach in the low 
mass range would be to exploit any future experimental facilities aiming to measure the 
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coherent neutrino-nucleus elastic cross-section. Indeed, it is well known that a source of 
stopped pions can be used to detect the elastic scattering of neutrinos on nuclei in DM-type 
detectors, sensitive to a recoil energy in the (9(10 keV) range. Should MeV DM exist, it would 
also produce a considerable recoil, possibly dominating over the coherent neutrino scattering 
signal. Such a possibility provides some additional physics motivation for the proposals such 



as CLEAR 41 , that are capable of detecting the coherent scattering of neutrinos. Finally, 
we note that modern long-baseline neutrino facilities such as MINOS and T2K have more 
energetic beams and may open the possibility of extending the mass reach at the upper end. 
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